A virtual intersubband spin-flip spin-orbit coupling induced spin relaxation in GaAs 

(110) quantum wells 
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A spin relaxation mechanism is proposed based on a second-order spin-flip intersubband spin-orbit 
coupling together with the spin-conserving scattering. The corresponding spin relaxation time is 
calculated via the Fermi golden rule. It is shown that this mechanism is important in symmetric 
GaAs (110) quantum wells with high impurity density. The dependences of the spin relaxation time 
on electron density, temperature and well width are studied with the underlying physics analyzed. 
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In recent years, semiconductor spintronics has be- 
come a focus of intense experimental and theoretical re- 
search P, i, H i, S, @j. One of the key factors for 
the design of the spin-based device is to understand the 
spin relaxation such that the information is well pre- 
served before required operations are completed. In n- 
type zinc-blende semiconductors, like GaAs, the lead- 
ing spin relaxation mechanism in most situation is the 
D'yakonov-PereP (DP) mechanism @, S, H, [HI, which 
is from the joint effects of the momentum scattering 
and the momentum-dependent effective magnetic field 
(inhomogenous broadening 

EL GJJI1) induced by the 
Dresselhaus [l5j] and/or Rashba 16. Il7| spin-orbit cou- 
pling (SOC). However, in symmetric (llO)-oriented GaAs 
quantum wells (QWs), when only the lowest subband is 
occupied, the in-plane components of the spin-orbit field 
vanish and the effective magnetic field only exists along 
the growth direction [li| . As a result, electrons with spin 
polarization along the growth direction can not precess 
around the effective field, and therefore the DP mecha- 
nism is absent. It is also noted that in the presence of 
an in-plane external magnetic field, the in-plane and out- 
of-plane spin relaxations are mixed, thus the DP mecha- 
nism still leads to the spin relaxation/dephasing in this 
system, as point out by Wu and Kuwata-Gonokami [lj| . 
Experimentally Ohno et al. [2fjj first observed very long 
spin relaxation time (SRT) in GaAs (110) QWs, which 
exceeds the SRT in (100) QWs by more than one order of 
magnitude. Later the spin dynamics in this system was 

Hi ni m m m m m, us m . 



studied by many works 



In most of these works, the main reason limiting the SRT 
is thought to be the Bir-Aronov-Pikus (BAP) mecha- 
nism [30l l3ll | , which is due to the electron-hole exchange 
interaction (holes are from the optical excitation in n- 
type samples). In the spin noise spectroscopy measure- 
ments by Miiller et al. 28], the excitation of semicon- 
ductor is negligible, and hence the BAP mechanism is 
avoided. They reported even longer SRT about 24 ns and 
attributed it to the DP mechanism due to the random 
Rashba fields caused by fluctuations in the donor den- 



sity 2S 33|. In this note, we propose another spin relax- 
ation mechanism to understand the slow spin relaxation 
in the absence of the DP and BAP mechanisms, which is 
based on a second-order spin-flip process of the intersub- 
band SOC together with the spin-conserving scattering. 
We will show that this mechanism is important in sym- 
metric GaAs (110) QWs with high impurity density. 

We start our investigation from symmetrically doped 
GaAs (110) QWs with the growth direction along the z- 
axis. The well width is assumed to be small enough so 
that only the lowest subband is occupied for the tem- 
perature and electron density we discuss. The envelope 
functions of the relevant subband are calculated under 
the finite-well-depth assumption. The barrier layer is 
chosen to be Alo.3gGao.6iAs where the barrier height is 
328 meV [34| . The Hamiltonian can be written as (H = 1) 



H = 



E 

nn'kerer 



[(e k + E n )5 n 



(k) 



(1) 



Here ek = k 2 /2m* is the energy spectrum of the elec- 
tron with two-dimensional momentum k = (k x ,k y ); E n 
represents the quantized energy of the electron in the n- 
th subband; <x is the Pauli matrix. In symmetric QWs 
without external gate voltage, the Rashba SOC can be 
neglected 35(, and the spin-obit field h(k) is only from 
the Dresselhaus term. In the (110) coordinate system, 
h(k) can be written as 



(k) = lD [-(kl + 2k 2 y ){n\k z \n') + (n\kl\n')] 
(k) = Ajuk x k y (n\k z \n'), 
h n z n '(k) = lu k x (kl-2k 2 y -(n\k 2 z \n>))6 nn ,, 



(2) 



where Tn = 30 eV-A 3 denotes the Dresselhaus SOC coef- 
ficient [36j] and {n\kf\n') = J dz^(z)(-id/d z ) m 4> n/ (z) 
with 4> n (z) representing the envelope function of 
the electron in rt-th subbands. Since (n\k z \n) = 
(n\k z \n) = 0, all the in-plane components of the in- 
tersubband spin-orbit field vanish as mentioned above. 
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However, (l\k z \2) ^ 0, thus the in-plane compo- 
nents of the intersubband spin-orbit field are still 
present. The interaction Hamiltonian Hj is composed 
of the electron-ionized-impurity, electron-longitudinal- 
optical (LO)-phonon, electron-acoustic (AC)-phonon and 
electron-electron Coulomb interactions. Their expres- 
sions can be found in textbooks 



37, 36 





FIG. 1: (Color online) Schematic representation of a second- 
order spin-flip process based on the intersubband spin-orbit 
field H so and the spin-conserving scattering H ac - 

A second-order spin-flip process can be constructed via 
the intersubband spin-orbit field and the spin-conserving 
scattering as depicted in Fig. 1. The matrix element of 
the spin-flip transition from the electron-impurity scat- 
tering (which is referred to as the spin-flip electron- 
impurity scattering in the following) is given by 



_ (lkt |g S o|2k|) (2k I \H s c f\lk> l) 

£k' - £k - AE 2 i 
(lktlTOk' T) (2k' t |ff„o|lk' |) 
-AE 21 



(3) 



Here AE21 = Ei — E\ stands for the energy splitting 
between the first and second subbands of electrons. Al- 
though the above process is similar to the short-range 
Elliott- Ycfct spin relaxation caused by virtual scatter- 
ing P, , the intermediate virtual states are chosen to 
be the states in the high conduction subband in stead of 
the states in the valence band in the previous work [3^ |. 
From the energy conservation e k = e k ', the symmetry 
of the form factor I\2{q z ) — hiilz) = (l\e lqzZ \2) and 
the symmetry of the spin-orbit field h 21 (k) = — h| 2 (k), 
C/jl(k, k', q z ) can be rewritten as 

U H (k,k', qz ) = ^||^ [h 12 (k)-h 12 (kO]-<r n . (4) 

Similarly, one can obtain the matrix element of the spin- 
flip electron-phonon scattering 

Dft (k, k',q z ) = D ^* { « z) [h 12 (k) - h 12 (k' )] • a n , (5) 

where Q = (k' — k, q z ) is the three-dimensional momen- 
tum, t/k'-k and Dq v represent the matrix elements of 
the electron-impurity and electron-phonon interactions 
respectively, whose expressions are given in detail in 



Ref. |4(J. From Eqs. (gj) and ©, it is seen that the 
matrix elements of the spin-flip electron-impurity and 
electron-phonon scatterings are proportional to |h 12 (k) — 
h 12 (k')|/A J B 2 i- Thus the term containing fcf in h(k), 
which is independent of two-dimensional momentum k, 
has no contribution to the spin-flip scattering. 

By using the Fermi golden rule, the SRT is given by 



ei/ep 



T ep- 



With \M 



2 £ r?f p (k,k', fe )/ k (i-/ k , 



kk'q z 



E /k(i - fw) 

kk' 



(6) 



Here / k is the equilibrium electron distribution func- 
tion. r» 1 / ep (k, k', q z ) represents the transition rate of the 
spin-flip electron-impurity or electron-phonon scattering, 
which can be written as 

rftCk.k'.g.) = 2^|{/ Ti (k,k',^)| 2 ^e k -e k 0, (7) 

r-J(k,k / ,fc) = 2^^|^ i (k,k', (?2 )| 2 [7V Q ^(e k -e k , 
n 

+ w Qr,) + (Nq v + l)5(e k - e k - - ujq v )] ■ (8) 

Here Nq^ and wq, are the distribution function and en- 
ergy spectrum of the phonon with mode rj and momen- 
tum Q. It is noted that the spin relaxation process we 
discuss is only from the virtual intersubband scattering. 
The spin relaxation due to the real intersubband scatter- 
ing 14 , 23, [25| is negligible due to the small well width 
in the present investigation. 

Our main results are summarized in Figs. 2-4. In 
Fig. 2, we plot the SRTs due to the spin-flip electron- 
impurity, electron-AC-phonon and electron-LO-phonon 
scatterings together with the total SRT with all the spin- 
flip scatterings included as function of temperature for 
Ni = N e (solid curves) and N, = 0.01 N e (dashed 
curves). N e = 1.8 x 10 11 cm -2 and a — 16.8 nm. 
In the low impurity density case with N% — 0.01 N e , 
it is seen that the spin relaxation due to the spin-flip 
electron-AC-phonon and spin-flip electron-impurity scat- 
terings are comparable at low temperature and the con- 
tribution from the spin-flip electron-LO-phonon scatter- 
ing becomes more important and even dominates at high 
temperature. In the high impurity density case with 
Ni = N e , the spin relaxation due to the spin-flip electron- 
impurity scattering is dominant at most temperatures we 
discuss, and the contribution from the spin-flip electron- 
LO-phonon scattering becomes comparable to that from 
the spin-flip electron-impurity scattering when T > 80 K. 

We also compare our results with the experiment by 
Mullcr et al. [28| (square) in Fig. 2. By fitting the mea- 
sured mobility given in Ref.l28l one obtains N ~ 0.01 N e . 
In this case, it is seen that the SRT obtained from our 
model is two orders of magnitude larger than the experi- 
mental data at T = 20 K. However, in the high impurity 
case with N = N e , it is shown that the spin relaxation 
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FIG. 2: (Color online) SRTs due to the spin-flip electron- 
impurity, electron- AC-phonon and electron-LO-phonon scat- 
terings as well as the total SRT with all the spin-flip scatter- 
ings included vs. temperature T with well width a = 16.8 nm 
and electron density N e = 1.8 x 10 11 cm -2 . The impurity 
densities are Ni = N e (solid curves) and 0.01 N e (dashed 
curves), respectively. The square represents the experimental 
result in Ref. [28]. The cross represents the SRT due to the 
RRDP mechanism for Ni = iV«, calculated via the relation 



oc 1/tj,. 
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FIG. 3: (Color online) SRTs with all the spin-flip scatterings 
included vs. temperature T for Ni = N e = 0.5 x 10 11 cm -2 
and 2 x 10 11 cm -2 , respectively. The well widths are a — 
10 nm (•) and 16.8 nm (□). The black dashed vertical lines 
indicate the Fermi temperatures T§ for both electron densi- 
ties. 



rate due to our mechanism is enhanced by more than one 
order of magnitude. Meanwhile, the spin relaxation rate 
due to the random Rashba field induced DP (RRDP) 
mechanism [32J, [33] is suppressed significantly due to the 
increase of the electron-impurity scattering. Here we es- 
timate the SRT due to the RRDP mechanism at Ni = N e 
(cross) by using the relation T^ RDF oc l/r p Q. It is seen 
that the mechanism from our model becomes more im- 



FIG. 4: (Color online) SRTs with all the spin-flip scatterings 
included vs. electron density N e for temperatures T = 20 K 
and 50 K, with impurity densities Ni = N e (•) and Ni — 
10 11 cm~ 2 (□). The black dashed vertical lines indicate the 
electron densities satisfying = feT for both temperatures. 



portant compared with the RRDP mechanism in the high 
impurity case. 

Now we discuss the temperature dependence of the 
SRT. In Fig. 3 we plot the SRT as function of temper- 
ature for different electron densities and well widths. It 
is shown that the SRT first decreases slowly and then 
rapidly with T. The turning point is roughly around 
the electron Fermi temperature Tp = Ep/ks- Since the 
electron-impurity scattering has a weak temperature de- 
pendence, the temperature dependence of the SRT is 
mainly from the fc-dependent spin-orbit field h(k). From 
Eqs. (|3]) and it is seen that only the quadratic terms 
in h(k) contribute to the spin-flip scattering. When tem- 
perature increases, more electrons are distributed at lager 
momentum states, and then the contribution from h(k) 
becomes larger. This leads to a larger spin-flip scatter- 
ing rate and hence reduces the SRT. It is also noted the 
average momentum of electrons is not sensitive to tem- 
perature in the degenerate regime (T < T|). Thus the 
SRT decreases slowly at low temperature. From Fig. 3, 
it is also seen that the SRT becomes shorter for wider 
well width. It can be understood as follows. Since the 
form factor is weakly dependent on well width, one ob- 
tains that T x oc (A£ 21 /(l|^|2)) 2 from Eqs. It 
is known that AE21 oc a" 2 and (l|fc 2 |2) oc a" 1 under the 
infinite well-depth assumption. Thus the SRT decreases 
with an increase of well width. 

The electron density dependence of the SRT is also in- 
vestigated. In Fig. 4, the SRT is plotted as function of 
electron density for different temperatures and impurity 
densities. We first discuss the case with fixed impurity 
density N — 10 11 cm~ 2 . It is seen that the SRT first in- 
creases and then decreases with iV e with a peak around 
the electron density satisfying Ep — k^T. This peak 
originates from the competition of two effects. On one 
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hand, similar to the discussion in temperature depen- 
dence, with the increase of electron density, the average 
k increases and thus the contribution from h(k) is en- 
hanced. This effect leads to a decrease of the SRT. On 
the other hand, the screening of electrons also increases 
with electron density, which suppresses the spin-flip scat- 
tering and increases the SRT. When the electron density 
is low and E§< < ks,T, i.e., in the nondegenerate regime, 
the average k changes little with N e . Thus the effect of 
the increase in the screening is dominant and the SRT 
increases with increasing density. However, when the 
electron density is high enough so that Ep > fceT, i.e., 
in the degenerate regime, the effect of the increase of the 
contribution from h(k) becomes more important. Conse- 
quently the SRT decreases. Then we turn to the case with 
Ni = N e . It is seen that the SRT decreases monotonically 
with electron density. This is because the effect of the in- 
crease of the impurity scattering surpasses the one of the 
increase of screening, and makes the SRT decrease even 
in low electron density (nondegenerate) regime. Never- 
theless, one still can see that the SRT decreases more 
rapidly in the degenerate regime. 

In conclusion, we have proposed a spin relaxation 
mechanism based on a second-order spin-flip process of 
the intersubband spin-flip SOC and the spin-conserving 
scattering, and calculated the corresponding SRT via the 
Fermi golden rule. We show that this mechanism is im- 
portant in symmetric GaAs (110) QWs with high impu- 
rity density. The temperature, well width and electron 
density dependences of the SRT are also investigated. 
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